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Thermodynamic Constraints on the Magnetic Field Dependence of the Neutron Resonance in 

Cuprate Superconductors 
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In order to test the recently proposed connection between the neutron resonance and the specific heat anomaly 
in cuprates, the experimental specific heat data on YBa.2Cu3O6.93 and a theoretical estimate of the single 
particle fermionic contribution to the specific heat is used to provide an upper bound for the intensity of the 
neutron peak. The deduced peak intensity is similar in magnitude and temperature dependence to that observed 
in neutron scattering experiments, and it is constrained to decrease strongly under the influence of moderate 
magnetic fields oriented along the c-axis of the crystal. An explanation is proposed for the predicted suppression, 
based on the observation that the resonance intensity is very sensitive to superconducting phase correlations. 



PACS numbers: 74.25.Bt,74.25.Ha,74.72.-h,74.20.Mn 

The resonance observed in inelastic neutron scattering re- 
mains one of the most heavily debated features of cuprate 
superconductivity. Its initial observation in YBa 2 CusOQ+ x 
[1] sparked numerous experimental [||] and theoretical 
1 3 — fT~r| I investigations. The detection of the resonance in 
Bt 2 Sr 2 CaCu 2 8+s together with several recent stud- 
ies suggesting a direct link between the resonance and salient 
features of other experiments [MOH], give further support to 
the suspicion, that the neutron resonance is a key phenomenon 
in high T c superconductors. 

Among the possible links, the one between the thermo- 
dynamic and magnetic anomalies deserves special attention. 
Taking a general argument by Scalapino and White fll5| ] fur- 
ther, Demler and Zhang [16| pointed out that the emergence 
of the neutron resonance below the superconducting transi- 
tion T c results in a decrease in exchange energy that could 
account for the entire superconducting condensation energy. 
They proposed a simple relation between the resonance in- 
tensity I Tes (T) and the corresponding resonance contribution 
C les (T) to the specific heat: 



Cres(T) 



J /dT res (T) 



2/4 



dT 



(1) 



Here /is is the Bohr magneton and J is the nearest-neighbor 
spin-spin coupling. Dai and his collaborators [|7|| calculated 
drlresiT) from neutron scattering data on Y Ba 2 CusOQ+ x , 
used Eq.(|l]) to compute the specific heat anomaly correspond- 
ing to the neutron resonance, and compared the outcome to 
the data of Loram and his collaborators [18]. They concluded 
with the remarkable observation that the data are compatible 
withEq.fg). 

However, if Eq.dTh holds, then its inverse must also be valid. 



Ires 

(T) = (2/4/J) / dtC Ie8 (t), 
Jt 



(2) 



where T is an upper cutoff temperature, chosen to be larger 
than T*, the resonance onset temperature. Thus, it should 
be possible to check the link between the thermodynamic 
and magnetic measurements by starting from the specific heat 



measurements instead. After only a qualitative examination 
of typical specific heat data C e i(T) (see, for example inset of 
Fig. [l]) one concludes that performing such a task is not possi- 
ble without an adequate account of the fermionic contribution 
to the specific heat. In fact, the specific heat jump in conven- 
tional superconductors is entirely due to the onset of a finite 
gap in the quasiparticle spectrum |[l9||. 

More important is the observation that if C res is defined ad- 
equately, the inverted relation Eq.(Q) allows one to impose a 
severe constraint on I res (T). It is known for some time, that 
the specific heat anomaly is removed [|20"|4?4"|] by a magnetic 
field of several Tesla oriented along the c-axis of the cuprate 
crystal. This means that Eq. (||), used with thermodynamic 
data taken in magnetic field, predicts large field induced ef- 
fects on the resonance intensity. This is a new property of the 
neutron resonance, that has not been observed yet. The aim 
of this paper is therefore to point out that (a) an upper bound 
can be given for C res (T) in Y Ba 2 Cu^O '5.93, which in turn 
yields a reasonable account of I reB (T), (b) the magnetic field 
sensitivity of C res (T) constrains I res (T) to be strongly field 
dependent, and (c) an intimate connection exists between the 
presence of a sharp resonance and superconducting long range 
order. This last observation leads to a natural explanation of 
the predicted field-induced suppression in resonance intensity. 

The procedure for obtaining an upper bound for C res (T) 
is based on the observation that one can give a lower bound 
for the specific heat contribution Cf(T) of single particle 
fermionic excitations. The fermionic specific heat contribu- 
tion can be calculated from fi/(T, V), corresponding to the 
fermionic part of thermodynamic potential [E5J] 



Vf(T,V) = --Tr 



— Sk^G 



<v«J> 



(3) 



where the single particle fermion propagator has the Nambu 
form Gk,f = (iQ — £k T 3 — ^kTi) 1 ^ corresponding to 
the Bardeen-Cooper-Schrieffer (BCS) self energy Sk,f = 
AkTim are Pauli matrices, Q = (21 + l)irk B T, and A k (T) 
is the superconducting gap. After performing the trace over 
the Pauli matrices, and the Matsubara sum, one obtains 
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n f (T,v) = {I Mi - M + ^(i - 2/0}, 

k,er 

(4) 

where /k = [1 + exp(/3i?k)]^ 1 is the Fermi function. This is, 
not surprisingly, the usual form of the BCS thermodynamic 
potential [^9|], since in a conventional superconductors the el- 
ementary excitation are fermionic in nature. 

The same statement does not necessarily hold for the 
cuprate superconductors. It was argued [[If]] that the specific 
heat jump in the cuprates is in large part due to the appearance 
of new bosonic degrees of freedom below T c , The degradation 
[ |l8| ] of the specific heat jump into an anomalous cusp in the 
underdoped regime adds further support to the possibility that 
jump is caused by the presence of bosonic excitation in the 
superconducting state. The expression in Eq. (Q) will be used 
to separate from C c \ the part that could, in principle, be due 
to bosonic excitations. This will be done by subtracting from 
Cci the fermionic contribution, computed from Eq. (Q) and the 
thermodynamic relation C* f (T) = -Td 2 n f (T, V)/{dT) 2 . 
There is no information available on the temperature depen- 
dence of Ak(T) in Y Ba-iCuzOs+x. Under these circum- 
stances it is more useful to provide a controlled lower bound 
for the fermionic specific heat. This is possible to achieve by 
taking the low temperature gap value A ~ 25 meV from tun- 
neling measurements [p6[], and calculating Cf (T) under the 
assumption that A is constant at all temperatures. In this 
way one obtains a lower bound for the fermionic specific heat 
contribution, as compared to the case when A(T) is decreas- 
ing with increasing temperature. Clearly the constant gap 
model systematically underestimates C/(T) in a system such 
as YBCL2CU3O6.93, where A(T) is known to close near T c . 
In the underdoped Y Ba^Cu-iO^x samples, due to the onset 
of the pseudogap, the constant gap model becomes an even 
better approximation [^7|| . Unfortunately, no tunneling data is 
available yet for such crystals. 

Within the constant gap model Eq. (Q) gives the following 
result for the fermionic specific heat 



Cf(T) 



fee 
~2~ 



^ sech 3 



,0E V 



(/3e k ) cosh( 



-(P 3 AiE k ) sinh( 
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The anomalous part, C rcs (T), will be defined as the differ- 
ence in the measured [ |T8| ] electronic specific heat C e \(T) and 
the fermionic contribution 



C ies (T) ~ Cei(T) - Cf (T). 



(6) 



The specific heat data of Loram et al. for YBa2Cu3C>Q.g3, 
displayed in in the inset of Fig.[lJ together with Eqs. ([|h(^]) 
-(||), can now be used to compute the temperature depen- 
dence of the neutron resonance, as predicted by the ther- 
modynamic data. At low temperatures this procedure gives 



-fres(O) ~ 0.01 /^f.u. -1 , to be compared to the correspond- 
ing experimental value of 0.06 ± 0.04/^f .u. _1 , where f.u. 
stands for formula unit. 

CJT) - C/T) [J/mole-K] 




FIG. 1. Inset:The measured specific heat C e i(T) (dots), and 
the lower bound for the fermionic contribution C/(T) (solid line), 
obtained by numerical evaluation of Eq.(^|), using an s-wave gap 
A = 25 meV at all temperatures. The normal state density of states 
was fitted to reproduce the normal state data near To. Main figure: 
The anomalous part, obtained by subtracting the calculated lower 
bound, Cf(T), from the experimental data, C e i(T) . The effect of 
the magnetic field along the c-axis of the crystal is simulated by using 
the results of Ref. [pl"|. The peak of the anomaly is gradually shifted 
to lower temperatures and reduced: 100% and 93 K (zero field, solid 
curve), to 71% and 89 K (B = 1 T, dotted curve), 45% and 81 K 
(B = 3T, dashed curve), and finally 30% and 75 K (B = 6 T, 
dash-dotted curve) 

While the deduced intensity is still the same order of mag- 
nitude as the experimental value of Ref. Jl7|], the agreement 
is less striking. On the other hand, the fact that the resonance 
intensity required by the thermodynamic changes is smaller 
than that observed in neutron scattering supports the possibil- 
ity that the source of superconducting condensation energy is 
in fact the gain in exchange energy, as suggested by Demler 
and Zhang [Jjl]] . The deduced I res (T) is somewhat affected by 
ambiguities in choosing To, by experimental errors in the tun- 
neling gap A, and in the specific heat data itself, introduced by 
the necessary phonon subtraction. The normalized zero field 
intensity is shown as the solid curve in Fig. ^, together with 
the experimental data from Ref. |17|]. The agreement is sat- 
isfactory except near the transition temperature T c ~ 92.5 K, 
where the intensity is overestimated. This can be considered 
an indirect signature of temperature dependence of the super- 
conducting gap. 

As pointed out by several groups J2C|-|2^], the specific heat 
anomaly near T c is removed in YBa2Cu30e.gs by a magnetic 
field B ~ 10 T along the c axis, while C e i(T) away from T c 
is unaffected. This is highly unusual, since the applied mag- 
netic field is only a small fraction of the upper critical field 
B C 2 ~ 100 T. Similar effects have been recently documented 
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in other cuprate materials as well [ |23|J24| ] . The recent work 
of Junod and coworkers [ E4| ] presents an especially detailed 
account of the field dependence of the specific heat anomaly. 
They contrast the behavior seen in the cuprates with that in 
a low T c superconductor (Nb-Zr alloy), where the anomaly 
is pushed to lower temperatures, but never removed. Indeed, 
the behavior seen in cuprates seems to go beyond the conven- 
tional framework, according to which as long as T c is finite, 
the specific heat jump is also finite [EJJ. 

The most detailed set of data for the field dependence of 
the anomaly in YBCI2CU3O6.93 was presented in Ref. pT|]. 
In these measurements, however, the phonon background was 
not subtracted, as in the data Loram et al. used above. In or- 
der to avoid subtraction problems, the finite field specific heat 
jumps were obtained by imposing on Loram and coworkers' 
zero field data [ |l8| ] the suppression in the height and the shift 
towards lower temperatures of the specific heat maximum ob- 
served by Junod's group J2l|]. The result is shown in the main 
part of Fig. |J 

When field-suppressed specific heat anomalies are used in 
Eq.(Q), the resulting temperature dependence of the deduced 
resonance intensity deviates from what was measured in zero 
field, as shown in Fig. |[ In fact, these results show that 
I 1CS (B,T) should be depleted at low temperatures beyond 
the present experimental error when B > 6 T (see figure). 
Bourges et al. [[l9| measured the neutron resonance in mag- 
netic field and found I 1CS (T ~ 6 K) to be insensitive to fields 
of order B ~ 12 T oriented along the ab plane of the sample. 
Their result does not invalidate the claim of Dai et al., since 
in that geometry the specific heat anomaly is only moderately 
depleted, and not removed by the field [Ell]. However, the 
measurement of I YOS (B, T) with B along the c-axis seems to 
be a crucial experiment for not only establishing the validity 
of Eqs.(JI[) and (||), but also to gain important insight into the 
nature of the neutron resonance itself. 

The strong magnetic field dependence of the resonance pre- 
dicted by the above procedure raises immediately a key ques- 
tion: What is the main mechanism responsible for inducing 
such a large change in the neutron intensity? In the remainder 
of this paper I will argue that there is an intimate connection 
between the resonance and correlations in the phase of the su- 
perconducting order parameter. This connection, besides of- 
fering a natural explanation for the field sensitivity of the res- 
onance, also seems to have direct relevance to some puzzles in 
the zero field properties of the resonance, like the existence of 
the resonance in bilayer materials only [[|], and the presence 
of a broad precursor resonance in the pseudogap state [|i7|]. 

There are three important clues to consider. First, irrespec- 
tive of doping, a sharp resonance is observed only below T c , 
after off diagonal long range order is established. Second, the 
magnetic field required to affect the resonance is predicted by 
the above analysis to be much smaller than B C 2- Scanning 
tunneling spectra [ 30 1 show that such fields do not broaden 
out the superconducting singularities significantly. Finally, 
a c-axis field strongly affects the specific heat jump and the 
resonance, but only mild changes are observed if the field is 



along the ab-plane [Elj], consistent with the results of Bourges 
et al. [EJJ. The strong dependence on the field direction sug- 
gests that Zeeman splitting is irrelevant and that the resonance 
is not affected by shielding currents along the c-axis itself, 
when the field is parallel to the CuOi planes. On the con- 
trary, the resonance is affected considerably by the presence 
of in-plane shielding currents circulating around the cores of 
Abrikosov vortices, when the field is parallel to the c-axis. As 
discussed in more detail below, this translates to the statement 
that in-plane long range coherence as well as short range c- 
axis coherence (intra-bilayer) does matter for the resonance, 
whereas long range c-axis phase coherence does not. 



O neutron scattering data 

From C,(T) data, B = T 

B = 1 T along c axis 

B = 3 T along c axis 

B = 6 T along c axis 
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FIG. 2. The neutron resonance intensity for the optimally doped 
compound YBa.2Cu3O6.93: Dots are the zero field results of Ref. 
Jiv[] ; lines were derived from specific heat data of Refs. [[H|^l| as 
explained in the text. All the curves were normalized with the zero 
field, low temperature value of the deduced intensity, O.Ol/igf .u. _1 . 

In order to make these observations more transparent, it is 
instructive to examine the bare spin susceptibility, where the 
phase dependence is indicated explicitly 



x°(Q,in n ) = x^ G (Q,in n ) 

p, itt n - iuj m )A(p, iu m ). 



p,m 



(7) 



Here Xgg(Q'*^«) anc ^ Xff(Q'*^«) are trie normai an d 
superconducting contributions to the bare susceptibility, 

X^ F (q,w„) = ^Ek/^k+qXCe +v n )Fk(Q), and il m = 
2mirkBT. A similar expression holds for x GG (Q, if2„), 
but with the Gor'kov propagators replaced by the diagonal 
part of the superconducting Green's function. The function 
A(p,ioj m ) is the Fourier and Matsubara transform of the 
charge displacement correlation function [p2[|, A(r — r', r — 



= ( T A< 



,i4>{r,T) p — ic/>(r' , 



)]), and (j)(r, t) is the phase of the 



superconducting order parameter at position r and imaginary 
time t. 

Equation (Q) can provide new insights even for the zero field 
case. In the absence of phase fluctuations, corresponding to 
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the BCS mean field level, A(q, iui m ) — /35(q)<5 TO) rj, and one 
obtains the standard result used as a starting point in essen- 
tially all theoretical descriptions of the resonance [pj-[TT{]. In 
zero external field the phase fluctuations are uniform in space 
and Gaussian around the mean phase <po, and consequently 
the correlator A(p, u>) is still dominated by its low frequency, 
long wavelength value A(0,0) = cxp[— (S<fi 2 )]. Thus, phase 
fluctuations reduce the contribution of the anomalous suscep- 
tibility XfF' an d consequently the resonance intensity J^-|TT|], 
since neither the bare, nor the RPA-dressed xgg can account 
on its own for the resonance. In single layer materials, due 
to the weak interplane coupling and strong two dimensional- 
ity, quantum phase fluctuations can be quite large, and con- 
sequently A very small. In contrast to single layer materials, 
the strong Josephson coupling between bilayers can stabilize 
the superconducting state with a rigid phase. This could be 
the main reason why the neutron resonance has only been ob- 
served in bilayer materials, and not in single layer materials, 
such as La2~ x Sr x Cu04 [0]. The broad maximum observed 
JT^ ] in the pseudogap regime of underdoped Y Ba^Cu^O^x 
could be the precursor of the resonance below T c , provided 
that strong pairing fluctuations are present in this regime J33|]. 

In the Abrikosov vortex phase induced by a field along the 
c-axis, the weight of A(p, to) is no longer concentrated at 
long wave lengths and low frequencies. As shown by Glaz- 
man and Koshelev [|34j], the vortex lattice destroys the off- 
diagonal long range order in the planes at finite T, while long 
range order along the field direction is broken only at ele- 
vated temperatures. For the in plane correlations, they find 
linv^oo A(r, z = 0) = lim r ^ 00 exp[—r/£(B)] = 0, where 
£(B) oc B^ 1 / 2 is the in-plane phase correlation length in 
magnetic field. Thus, the field will induce a momentum broad- 
ening <5p oc \J~B in the resonance, and a reduced intensity. 
In general, a form of A(q, lo) that deviates from the singular 
mean field result, will upset the resonance condition satisfied 
by the zero field susceptibility for Q = (ir, ir, ir); lo — loq. 

The arguments presented above are quite general and con- 
sequently have direct relevance not only to the magnetic field 
dependence of the neutron resonance, but also to the finite 
field behavior of optical conductivity, tunneling density of 
states, Raman scattering cross section, and other spectro- 
scopic and thermodynamic properties of cuprate superconduc- 
tors. A detailed study of all these effects is currently under 
way, and will be reported in a later publication. 

I am grateful to Dr. John W. Loram for making his spe- 
cific heat data available. I would also like to thank Prof. A. 
A. Abrikosov, Dr. Joel Mesot, and especially Dr. Michael 
R. Norman for numerous discussions, and to Dr. Philippe 
Bourges, Dr. Pengcheng Dai and Dr. Herbert A. Mook for 
correspondence. This research was supported in part by the 
NSF under awards DMR9 1 -20000 (administrated through the 
Science and Technology Center for Superconductivity), and 
the U.S. DOE, BES, under Contract No. W-31-109-ENG-38. 



[9 
[10 
[11 
[12 
[13 

[14 



[15 
[16 
[17 
[18 
[19 

[20 
[21 

[22 

[23 
[24 
[25 



[26 

[27 
[28 



[29 
[30 



[31 

[32 
[33 
[34 



J. Rossat-Mignod et al. Physica C, 185-189, 86 (1991). 

For a recent review, see T.E. Mason in Handbook on the Physics 

and Chemistry of Rare Earth s, Edited by K. A. Gschneider Jr. et 

al. dcond-mat/ 98 122 871). 



M.Lavagna and G. Stemmann, Phys. Rev. B 49, 4235 (1994) 
S. Quinlan and D.J. Scalapino, Phys. Rev. B 51, 497 (1995); N. 
Bulut and D.J. Scalapino, Phys. Rev. B 51, 5149 (1996) 
E. Demler and S.-C. Zhang, Phys. Rev. Lett. 75, 4126 (1995); 
Phys. Rev. B 58, 5719 (1998). 

D.Z. Liu, Y. Zha and K. Levin, Phys. Rev. Lett. 75, 4130 (1995). 
I.I. Mazin and V. Yakovenko, Phys. Rev. Lett. 75, 4134 (1995). 
L. Lin, S. Chakravarty, and P.W. Anderson, Phys. Rev. Lett. 78, 
3559 (1997). 

A.A. Abrikosov, Phys. Rev. B 57, 8656 (1998). 
D.K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 (1998). 
J. Brinckmann and P.A. Lee, Phys. Rev. Lett. 82, 2915 (1999). 
H. Fong et al., Nature 398, 588 (1999). 

J.C. Campuzano et al, Phys. Rev. Lett. 83, 3709 (1999); Ar. 
Abanov and A.V. Chubukov, Phys. Rev. Lett. 83, 1652 (1999). 
D. Munzar, C. Bernhard, and M. Cardona, Physica C 318, 547 
(1999). J.P. Carbotte, E. Schachinger and D.N. Basov, Nature 
401, 354 (1999). 

D. J. Scalapino and S.R. White, Phys. Rev. B 58, 1347 (1998). 

E. Demler and S.-C. Zhang, Nature 396, 733 (1999). 
P. Dai et al, Science 284, 1344 (1999). 

J.W. Loram et al, J. Supercond. 7, 243 (1994). 

G. Rickayzen, Theory of Superconductivity, John Wiley & 
Sons, NY, 1964. 

S.E. Inderhees et al, Phys. Rev. Lett. 66, 232 (1991). 

E. Janod et al, Physica C, 234, 269 (1994). M. Roulin, A. 

Junod, and E. Walker, Physica C 260, 257 (1996). 

N. Overend, M.A. Howson, I.D. Lawrie, Phys. Rev. Lett. 72, 

3238 (1994). 

A. Carrington et al, Phys. Rev. B 55, R8674 (1997). 
A. Junod, A. Erb and Ch. Renner, Physica C 318 333 (1999). 
A.A. Abrikosov, L.P Gor'kov and I.E. Dzyaloshinski, Methods 
of Quantum Field Theory in Statistical Physics, Prentice Hall, 
Englewood Cliffs, N.J. 1963. 

J.M. Valles et al, Phys. Rev. B 44, 11986 (1994); J. Y. T. Wei, 
et al, Phys. Rev. Lett. 81, 2542 (1998). 

H. Ding et al, Nature 382, 51 (1996). 

A. L. Fetter and PC. Hohenberg, Theory of Type 11 supercon- 
ductors, Vol. II, Chapter 14 of Superconductivity Edited by R.D. 
Parks (Dekker, NY 1969). 

P. Bourges et al, Physica B, 234-236, 830 (1997). 
Ch. Renner et al, Phys. Rev. Lett. 80, 149 (1998); E. D. Hudson 
etal, Science 285, 88 (1999); A. Yazdani etal, Phys. Rev. Lett. 
83, 176 (1999). 

The corresponding expression for the c-axis conductivity has 
recently been suggested by L.B. Ioffe and A.J. Millis, Science 
285, 1241 (1999). 

See, for example B. Janko and V. Ambegaokar, Phys. Rev. Lett. 
75, 1154 (1995) and references therein. 

See J. Maly, B. Janko and K. Levin, Physica C 321, 113 (1999), 
and references therein. 

L.I. Glazman, and A.E. Koshelev, Phys. Rev. B 43, 2835 (1991). 



4 



